SU(3) glueballs on coarse, anisotropic lattices* 
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Results from a calculation of the low-lying glueball spectrum of pure-gauge SU(3) are used as a test of the 
effectiveness of improved discretisation schemes in reducing finite spacing errors. Glueball masses are extracted 
from simulations on anisotropic lattices, where the temporal lattice spacing is much shorter than the spatial ones. 
This allows clearer resolution of the decay of glueball correlators. 



1. INTRODUCTION 

Current interest in the use of Symanzik's im- 
provement scheme, supplemented by mean-field 
link renormalisation, to provide a better discreti- 
sation scheme for QCD has stemmed from recent 
successes of coarse lattice simulations jlj . Coarse 
lattices offer the advantage of a significant reduc- 
tion in computational overheads, since the ap- 
proach to the continuum limit is accompanied 
by a rapid increase in CPU time requirements 
for Monte Carlo calculations. The improvement 
scheme allowed, for example, calculations of the 
static potential at lattice spacings up to w 0.4 fm 
with discretisation errors of only 5%. An attempt 
to use these actions on coarse lattices to exam- 
ine the glueball spectrum of QCD proved diffi- 
cult, however [§. The QCD glueballs arc mas- 
sive states and their evaluation on the lattice is 
notoriously difficult due to the large vacuum fluc- 
tuations of glueball operators. On a lattice with 
coarse temporal discretisation, the glueball cor- 
relator can barely be resolved before it is lost in 
the noisy vacuum. The use of a large operator ba- 
sis set to optimise the ground state overlap was 
hampered by a technical difficulty caused by the 
2x1 rectangle in the improved action. This term 
couples next-to-nearest neighbouring time-slices 
and thus modifies the transfer matrix, making 
it non-hermitian (physically, this term leads to 
the presence of unwanted extra modes in the im- 
proved gluon dispersion relation). A hcrmitian, 
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positive transfer matrix is required for the validity 
of the variational calculation employed and to en- 
sure that the effective mass approaches its asymp- 
totic value from above. These problems naturally 
suggest breaking the explicit Euclidean symmetry 
of the action and using a different lattice spac- 
ing in the spatial and temporal directions. With 
a shorter temporal lattice spacing, a t , clear evi- 
dence for plateaux in the effective masses can be 
seen before the signal is lost. Once discretisations 
of this form are used, a solution to the transfer 
matrix problem arises too. 

In this article, a study of three of the lighter 
glueball states at six different lattice spacings is 
discussed. For all but the scalar glueball, the re- 
sults show the expected scaling behaviour. 

2. ANISOTROPIC LATTICE QCD 

The improvement scheme can be extended 
readily to anisotropic lattices In this study, 
the following action was used: 
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with P„„ the plaquette in the fi— v plane, and R^ u 
the 2x1/1-;/ plane rectangle (with longest side 
parallel to the [i axis). This action, intended for 
use with a t <C a s , has 0(af, of) discretisation er- 
rors. The coefficients were determined using tree- 



level perturbation theory and tadpole improve- 
ment [0]. The spatial mean link, u s , was tuned 
such that the input value in the action matched 
its measured value as defined by u s — {Pss') 1 ^ 4 - 
This tuning is rapidly convergent and requires a 
minimal amount of work. The temporal mean 
link, u t , was fixed to Ut = 1, since its value in 
Landau gauge differs from unity by 0{a 2 /a 2 s ). 

3. SIMULATION DETAILS 

A set of six simulations were performed at fixed 
anisotropy, a t /a s = 1/3, with the intention of 
examining the scaling behaviour of the glueball 
masses using this action. The run parameters 
are given in Table |l|. Configurations were gener- 
ated using the Cabibbo-Marinari (CM) pseudo- 
heatbath and Creutz over-relaxation (OR) meth- 
ods. In our previous coarse lattice calculation g], 
large statistical samples were gathered in an at- 
tempt to resolve the glueball correlators. For this 
work, much smaller statistical samples were re- 
quired and the largest data set collected (for the 
coarsest lattice) consisted of 100 bins, with each 
bin containing an average of 100 measurements, 
each measurement separated by 3 sweeps of a 1-3 
hybrid CM/OR update. This is to be contrasted 
with the 334 bins of 1000 measurements from the 
old calculation, which yielded far less informa- 
tion. 

3.1. Setting the scale 

The lattice spacing was set by fitting the on- 
axis static potential to a standard ansatz, 

V(r) =V --+ar, (2) 
r 

then using the fit to determine the hadronic scale, 
ro 0, defined as r 2 dV/dr — 1.65 at r — ro, 
which corresponds roughly to ro ~ 0.5fm. This 
definition gives a lattice scale that is less depen- 
dent on the choice of potential fit ansatz than the 
string tension scale. The results of this determi- 
nation are given in Table |l|. In all these calcula- 
tions, the input ratio a t /a s has been used since 
its renormalisation is known to be small ||. In 
all runs, the lattice volumes were chosen to be 
large enough such that the masses of the scalar 
and tensor states should lie within i% of their 



Table 1 

Run parameters for our simulations, a t /a s = 1/3. 
The lattice spacing is set using ro = 0.53 fm. 
Errors in a s are purely statistical. 
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Lattice 


U s 


a s (fm) 


L(fm) 


1.7 


6 3 x 18 


0.745 


0.457(4) 


2.7 


1.9 


6 3 x 18 


0.764 


0.414(3) 


2.5 


2.0 


8 3 x 24 


0.772 


0.387(3) 


3.1 


2.2 


8 3 x 24 


0.789 


0.329(2) 


2.6 


2.4 


8 3 x 24 


0.806 


0.272(4) 


2.2 


2.6 


10 3 x 30 


0.819 


0.217(3) 


2.2 



infinite volume limits as determined by Luscher's 
finite volume analysis p]. 

3.2. Optimising the ground state overlap 

The low signal-to-noise ratio of glueball calcu- 
lations makes it crucial to find lattice operators 
with good ground state overlaps. Large bases 
of at least 24 operators were used and optimal 
correlators were determined using the variational 
method. The glueball operators were built from 
spatial Wilson loops constructed from either APE 
smeared or Teper fuzzed links 0. As in previ- 
ous investigations, the smearing and fuzzing tech- 
niques were found to be far more efficient when 
they included a gauge-invariant projection of the 
new links back onto SU(3). In most cases, the 
APE smeared operators gave larger contributions 
to the ground state; the glueballs extend over a 
small number of sites of the lattice, and thus oper- 
ators built from fuzzed "superlinks" are too large. 
At most, one iteration of Teper fuzzing proved 
useful. In all cases, these procedures gave excel- 
lent ground state overlaps and in many cases, the 
overlaps were statistically consistent with unity. 

4. RESULTS 

4.1. The scalar glueball 

The effective mass for the scalar glueball for 
one of the simulations performed is shown in Fig- 
ure 0. For this example, a correlator signal can be 
seen out to time-slice 7, and a good single expo- 
nential fit in the range 1 — 7 can be performed. A 
signal for the first excited state can be obtained 
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Figure 1. The effective mass for the scalar glue- 
ball at = 2.4 (a s » 0.27 fm). The lines show 
results from fitting the correlator to a single ex- 
ponential. 
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Figure 2. The effective masses for the tensor and 
pseudovector states at (3 = 2.4. The three irreps 
£ , ++,r 2 ++ and 7\ + ~ are plotted as □,<> and *, 
respectively. Solid lines show results from fitting 
correlators to single exponentials. 



from the first excited variational eigenvector. The 
mass of this state is slightly less than twice that 
of the ground state, consistent with an excited 
state observed in Wilson action calculations ||. 
The dependence on the lattice spacing is shown 
in Figure |^. The scalar mass shows a significant 
"dip" as the lattice spacing is increased, at the 
bottom of which the glueball is about 70% of its 
continuum value (from Wilson data). This dip 
is caused possibly by the presence of a critical 
point in the fundamental/adjoint plane, similar 
to that which affects the scalar glueball mass cal- 
culated using the Wilson action ||. No contin- 
uum extrapolations to the scalar data have been 
included. 

4.2. Other glueball states 

Simulation data for the tensor (2 ++ ) and pseu- 
dovector (1 H ) states were also examined. The 
(3 = 2.4 effective masses for the three lattice irreps 
(the T 2 ++ and E++ for the tensor and the for 
the pseudovector) with these continuum quantum 
numbers are shown in Figure [| Again, reliable 
single exponential fits to these correlators can be 



performed over a range of time-slices. The scaling 
behaviours of these states are included in Figure 
H For the E ++ and T± channels, no scaling vi- 
olations are seen while the T++ mass does show 
some cut-off dependence. At the largest lattice 
spacing, this leads to a 15% split in the masses, 
caused by rotational invariance breaking. The 
cut-off dependence is consistent with the antic- 
ipated leading discretisation error, proportional 
to (however other functional forms can not be 
ruled out). Extrapolations to the continuum limit 
are summarised in Table ||. Here, the simplest 
functional form that gives a reasonable x 2 is used 
in the extrapolations. Note that in the continuum 
limit, the masses of the two tensor irreps are in 
excellent agreement, suggesting restoration of ro- 
tational symmetry. 

5. CONCLUSIONS 

The simulations discussed here clearly demon- 
strate the advantages of using anisotropic actions 
to study the glueball spectrum of QCD. Mean 
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Figure 3. Scaling behaviour of the glueball states. 
The lattice irreps E++ , T++ and T+~ are 

labelled o, □, O and *, respectively. Crosses indi- 
cate Wilson action data from || . The solid line is 
a fit to the scalar glueball Wilson data, the dashed 
lines are fits to the leading scaling behaviour from 
the improved action (see Table ||) . 



link improvement was crucial in setting the cou- 
plings in the action. Using only workstations and 
an improved action with anisotropy 1/3, clear 
plateaux for the scalar, tensor and pseudovec- 
tor glueball states have been resolved, a pos- 
sibility previously reserved only for large-scale, 
supercomputer-based calculations. Including all 
six simulations, a total of approximately 5 x 10 9 
link updates were performed, two to three orders 
of magnitude fewer than Wilson action calcula- 
tions of similar statistical accuracy. 

Except for the scalar glueball, all the states 
demonstrate the expected scaling behaviour for 
this action. Using l/r = 372 MeV to set 
the scale gives a tensor mass of 2140 ± 45 MeV 
(from the E ++ irrep) and a pseudovector mass 
of 2620 ± 60 MeV. The errors quoted are purely 
statistical. In the case of the scalar, the improve- 
ment has reduced the cutoff contamination signif- 
icantly. 



Table 2 

Tensor and pseudovector continuum extrapola- 
tions. The data are fit to one of two alternatives: 
(A) m g (a s )r Q = c, or (B) m g (a s )r Q = c + da*. 





Fit type 


X Vdof. 


m g ro (continuum) 


E++ 


A 


0.52 


5.74 ±0.12 


rp + + 


B 


0.38 


5.77 ±0.08 




A 


0.29 


7.04 ±0.17 



Alternative improved actions which could have 
less contamination from the adjoint/fundamental 
fixed point are under investigation. With the suc- 
cess in extracting the continuum masses of the 
three states considered, an extended calculation 
to compute the masses of the lightest states in all 
20 lattice irreps seems feasible, perhaps using a 
higher anisotropy. 
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